Magnetic properties of the coupled ladder system MgV 2 Os 
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We present magnetic susceptibility measurements on MgV^Os, a compound in which the vanadium 
oxide planes have the same topology as in CaV2 0s. The most striking property is that there is 
an energy gap of about 15 K, much smaller than in CaX^Os, where the values reported are of the 
order of 500 K. We show that this may be understood naturally in terms of the phase diagram of 
the frustrated, coupled ladder system. This analysis leads to the prediction that MgV2 0s should 
have incommensurate dynamic spin fluctuations. 

PACS Nos : 75.10.jm 75.40.Cx 75.50.Ee 



The study of low dimensional, spin- 1/2 magnets is a 
very active field of research. One of the recent successes 
has been the synthesis and theoretical understanding of 
spin ladder systems jij . These consist of weakly interact- 
ing pairs of antiferromagnetic chains coupled by antifer- 
romagnetic rungs, and the single most interesting result 
is the appearance of a gap in the spin excitation spec- 
trum. The prototype of such systems is SrCu203, whose 
magnetic properties can be accounted for by a spin-1/2 
Heisenberg model 
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restricted to the exchange integrals depicted in Fig. 1. In 
the case of SrCu203, it is believed that | Ji| <C J2, J2, so 
that the system is clearly in the ladder limit. The fam- 
ily of planar vanadates also contains compounds where 
V atoms in ionization state V + give rise to spin-1/2 
planes having the topology of Fig. 1, and their primary 
difference from SrCu203 is that the coupling between the 
ladders is not a priori a small parameter. The first ex- 
ample in which a spin gap was reported is CaV205 
The presence of a large spin gap (about 500 K) in spite of 
non-negligible inter-ladder coupling came as a surprise, 
and prompted four of the present authors to study the 
phase diagram of the model depicted in Fig. 1 H. In 
particular, they showed that the ladders behave in an 
essentially decoupled manner when the second-neighbor 
coupling constants J2 and J2 are large enough compared 
to Ji, and that there should be an ordered, incommen- 
surate phase when the three coupling constants are com- 
parable. These results provide a reasonable explanation 
of the presence of a large spin gap in CaV205 since, as 
we will argue below, the information currently available 
on the magnetic interaction parameters leads to the es- 
timates J2/J1 — 3.5 and J 2 /Ji — 4.3 for that system. 



Materials which could be used to probe other regions of 
the rather rich phase diagram of Ref. || would be very 
welcome. 

In this paper, we present magnetic susceptibility mea- 
surements for another member of the vanadate family 
with the same topology, namely MgV2C>5. This system 
is quite intriguing because, although the vanadium ox- 
ide planes are very similar to those in CaV2C>5 [0, the 
magnetic properties are strikingly different, with, in par- 
ticular, a spin gap of only 15 K. As we shall argue below, 
this difference may be attributed to the presence of a he- 
lical ordered phase in the phase diagram of the model of 
Fig. 1, thus lending further support to the results of Ref. 

§• 

MgV2C>5 was prepared by solid state reaction, start- 
ing from a mixture of MgO and VO2, according to the 
procedure described in Ref. , and its structure was de- 
termined from x-ray powder diffraction Rietveld anal- 
ysis |J. It crystallizes in the orthorhombic system of 
space group Cmcm with, at 294 K, a = 3.6913(2)A, b = 
9.9710(4)A, and c = 11.0173(4)1. The structure retains 
the basic framework of the V2O5 host lattice, which is 
built up in one direction by infinite double lines of VO5 
square pyramids joined by alternate sharing of edges and 
corners to form vanadium zigzag chains along the a axis, 
and held together by corner sharing in the [001] direc- 
tion. It is therefore indeed analogous to the structure 
of the compounds a'-NaV 2 5 @ and CaV 2 5 [§. A 
peculiar feature of the MgV20s structure is the shift 
by a/2 along [100] for alternate layers (doubling of the 
short inter-layer parameter), and the puckering of the 
V2O5 layers (angle \x = 21°) to accommodate Mg ions in 
tetrahedral coordination. An idealized view of the struc- 
ture viewed down the a axis is presented in Fig. 2. As 
for CaV205, MgA^Os is a monovalent (V 4+ ) oxide, and 
vanadium atoms are located on a single crystallographic 
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site. 

Magnetic susceptibility measurements were performed 
using a superconducting quantum interference device 
QUANTUM DESIGN SQUID susceptometer. The mag- 
netic field induction was lkG. The molar susceptibility 
was corrected for diamagnetism by using Pascal's con- 
stants. The raw data are presented in Fig. 3. These data 
were reproduced on three different samples, with, how- 
ever, a large variation of the temperature-independent 
background from sample to sample. The main features 
of these data are: i) a regular increase between 400 K and 
100 K; ii) a small, not fully reproducible bump around 
85 K; iii) a very slow decrease between 70 K and 15 K; 
iv) a much faster decrease down to 4 K; v) a very fast 
increase at lower temperature. 

The low-temperature data are typical of a system with 
a spin gap. We have thus performed a fit of the data 
between 2 K and 25 K with the formula 

x(T) = A + B/T+(C/VT)exp(-A/T), (2) 

where A includes all temperature-independent contribu- 
tions (Van Vleck, plus possible contributions from impu- 
rity phases with temperature-independent susceptibility 
in this range), B/T accounts for paramagnetic impuri- 
ties, and (C/VT) exp(— A/T) describes the susceptibil- 
ity of a one-dimensional system with a gap || . The fit is 
depicted in the inset of Fig. 3, and the parameters used 
are A = 1.71 10 -3 emu/mole, B = 1.16 10 -3 emuK/molc, 
C = 7.2810- 3 emuK 1 / 2 /mole and A = 14.8AT. The fit 
is very accurate, with a residue of only 1.6 x 10 -6 over 
a temperature range 2 - 25 K. We have also tried to fit 
the data with a term of the form Cexp(— A/T) typical 
of gapped, two-dimensional systems J8| . The result was 
clearly not as good, which we take to indicate the effec- 
tive one-dimensionality of the system, or that MgV2C>5 
remains in the limit of weakly interacting ladders. 

While further investigation of this system, for exam- 
ple by Nuclear Magnetic Resonance, would be very use- 
ful to confirm these results, the presence of a rather 
small gap (~ 15K) as compared to CaV2C>5 (~ 50QK) 
is clearly very interesting. We believe that this result 
can be understood in terms of the model of Fig. 1, due 
to modest but not negligible differences in the structural 
parameters leading to important differences in the ex- 
change integrals. The direct determination of exchange 
integrals for a particular system is a very difficult task, 
because it requires an accurate fit of the temperature- 
dependence of the susceptibility over a large tempera- 
ture range. In vanadates this has been possible only in 
two cases for which we have a good theoretical under- 
standing of the temperature-dependence of the suscep- 
tibility of the underlying model, a'-NaX^Os pJlCl and 
CsV 2 5 Q. a'-NaV 2 5 is a one-dimensional Heisen- 
berg chain of corner-sharing VO5 pyramids with a single 
exchange integral J2(NaV2C>5). The susceptibility of this 
model is known very accurately from the work of many 
authors and a fit of the experimental susceptibility 



gives J 2 (NaV20s)~ 530K. CSV2O5 consists of weakly 
interacting pairs of edge-sharing pyramids. The relevant 
model is a Heisenberg model with only two spins coupled 
by an exchange integral Ji(CsV205). The susceptibility 
can be calculated analytically, and a fit to the experi- 
mental results gave Ji(CsV 2 5 )~ U6K 0. 

To obtain estimates of the exchange integrals in other 
members of the family, we make the reasonable assump- 
tion that these depend only on the local geometry of the 
bonds. For corner-sharing pyramids ^2^2), the dom- 
inant process is superexchange, and the exchange inte- 
grals are expected to depend on the local parameters de- 
fined in Fig. 4 according to || 

J 2 (d,9), J' 2 {d,9) cx cos 4 ((9)cr 14 , (3) 

where the exponent is given by empirical laws concerning 
the dependence of the overlap integrals on distance [jl3| . 
For edge-sharing pyramids (Ji), the origin of the inter- 
action is not so clear. Assuming that it is due primarily 
to direct exchange between vanadium 3d orbitals, similar 
empirical laws predict that it should scale according to 

Mdv-vj^idv-vr 10 - (4) 

where dy-v is the vanadium-vanadium distance. Having 
no other, more direct source of information at hand for 
vanadates, we will use these dependences in the following. 
The results for the exchange integrals of CaV2C>5 and 
MgV2C>5 based on the most accurate structural informa- 
tion available for a'-NaV 2 5 @, CsV 2 5 0, CaV 2 5 
H and MgV 2 5 @ are given in Table I. The trend is 
quite clear: Ji is larger for MgV205 than for CaX^Os, 
while J2 is considerably smaller, so that the ratio J2/J1 is 
approximately 4.3 for CaV20s, but only 2.5 for MgV205. 
According to the analysis of the model performed in Ref. 
||, this ratio is a very important parameter: J2/J1 = 
4.3 clearly puts CaV2C>5 in the ladder limit with a large 
gap |iq] . However, J 2 /Ji = 2.5 with J2/J1 = 2.8 puts 
MgV2C>5 very close to the helical ordered - hence gap- 
less - phase according to the phase diagram obtained by 
Schwinger boson mean-field theory The system is 
thus expected to have a very small gap, as this vanishes 
continuously at the boundary. Of course, estimates based 
on Eqs. (3) and (4) are not very accurate ]l7j| , but the 
conclusions can be expected to remain qualitatively cor- 
rect as long as the exchange integrals depend strongly 
on the distances, which has indeed been found to be the 
case in insulating oxides JlS| ]. 

An interesting consequence of the proximity to a heli- 
cal ordered phase is that MgV2C>5 should exhibit incom- 
mensurate dynamic spin fluctuations. We do not present 
calculations to elaborate on this point here. This pre- 
diction could be verified by inelastic neutron scattering 
when single crystals of sufficient sizes become available. 

Let us now turn to the high-temperature data. The 
increase of susceptibility with decreasing temperature, 
showing a maximum at intermediate temperatures, is 
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typical of systems with antiferromagnetic exchanges. Fit- 
ting the measured curve with a Heisenberg model has 
proven to be as difficult here as for other members of the 
family. The problem of the peak occurring well below the 
energy scale of the interactions may be due to poor theo- 
retical estimates of the temperature dependence for frus- 
trated systems with dynamic helical fluctuations. Work 
is in progress in this direction. 

We focus finally on the broad susceptibility peak which 
occurs around 85 K. To address the possibility that this 
behavior is related to structural modifications, we have 
studied the thermal evolution of the x-ray diffractograms. 
X-ray powder patterns were collected up to 8 — 30° (ex- 
cept the diffractogram at 83 K which was collected up 
to 75° in order to perform a Rietveld refinement |J) in 
a high-accuracy MICROCONTROLE diffractometer, us- 
ing CuKa radiation (graphite monochromator) from an 
18kW rotating- anode generator and a flowing cryostat 
with a stability range of 0.5 K. The thermal variation 
of the cell parameters, determined by least-squares re- 
finement, is presented in Fig. 5. The lattice parameter 
a along the vanadium chains increases slightly from 170 
K to 100 K, a temperature at which the fluctuations in 
the data suggest that a structural instability may be oc- 
curring, then decreases monotonically from 65 K to 8 K. 
Interestingly, the instability temperature range 100-65 K 
corresponds to the maximum observed in the susceptibil- 
ity curve. Moreover, a shortening of the lattice parameter 
a is qualitatively consistent with an increase of the ex- 
change integrals J\ and J2, hence with the experimental 
decrease of the susceptibility. However, it is difficult to 
correlate this with a structural modification because we 
did not observe a notable change on the Rietveld refine- 
ment of the structure performed at 83 K ||. The in- 
terlayer distance b decreases continuously over the whole 
temperature range, while c presents no clear evolution. 
A structural determination at low temperatures (around 
4 K) would be very useful, and we plan to perform this 
in the near future by neutron powder diffraction. 

To summarize, the magnetic susceptibility data re- 
ported in this paper suggest that MgV2 05 is a new spin 
gap system with a gap of about 15 K. This very small 
value should be contrasted with the much larger value 
(~500 K) reported for CaV2 05. We have shown that 
this difference can be understood in terms of the frus- 
trated, coupled ladder system due to the proximity to an 
ordered helical phase in the case of MgV20s. This in- 
terpretation suggests that MgV20s should have incom- 
mensurate dynamic spin fluctuations, a prediction which 
could be tested by inelastic neutron scattering as soon as 
suitably large crystals are available. Given the potential 
interest of systems with spin gaps in the range 0-30 K 
for studies under magnetic field, the magnetic properties 
of this system are likely to become a very active field of 
research. 
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3.0257 
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730 


3.5 


4.3 


MgV 2 5 


2.976 


201 


0.666 


1.954 


1.975 


565 


511 


2.8 


2.5 



TABLE I. Structural parameters and measured or esti- 
mated exchange integrals of CSV2O5, a'-NaV^Os, CaV^Os 
and MgVaOs. Distances are in A, and exchange integrals are 
in Kelvin. 
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FIG. 1. Schematic representation of the frustrated coupled 
ladder system with one type of nearest neighbour exchange in- 
tegral (Ji) and two types of next nearest neighbour exchange 
integrals {J21J2) 

FIG. 2. Idealized structure of MgX^Os viewed down the a 
axis (3.69A). 
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FIG. 3. Thermal variation of the molar magnetic suscep- 
tibility of MgV205. Diamonds: experimental points. Inset: 
enlargement of the low-temperature region with fit (solid line) 
to Eq. (2). 

FIG. 4. a) V-3d and 0-2p orbitals involved in superex- 
change between corner-sharing VO5 square pyramids; b) basic 
parameters of the VO5 square pyramid. 

FIG. 5. Temperature dependence of lattice constants for 
MgV^Os. Triangles and black dots correspond to two suc- 
cessive collection sets. The open circles correspond to the 
Rietveld refinement at 83 K. 



